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Final  Technical  Report  for  AFOSR  contract  89NM605 

During  the  reporting  period  (94/95)  very  good  progress  has  been  made  on  the  projects 
covered  by  this  contract. 

Collaboration  has  been  continued  with  Kurt  Bryan  related  to  crack  determination. 
Following  up  on  our  earlier  results  concerning  the  determination  of  multiple  cracks  from 
boundary  measurements  we  have  completed  a  study  of  the  relation  between  the  electrical 
boundary  data  coming  from  a  cluster  of  small  cracks  and  that  coming  from  a  single  effective 
macrosopic  crack  [1].  As  a  simple  example,  consider  a  two  dimensional  array  of  cracks  lying 
one  after  the  other  on  a  family  of  closely  spaced  parallel  lines.  In  this  case  our  analysis 
shows  that  if  the  “in  line”  distance  between  successive  cracks  is  exponentially  small  then 
the  data  looks  like  that  coming  from  a  single  macroscopic  crack.  If  not,  then  the  data  may 
be  used  to  distinguish  the  cracks!  In  addition  to  being  an  interesting  “homogenization”  - 
result,  this  also  provides  an  immediate  explanation  why  closely  spaced  “in  line”  cracks  are 
quite  accurately  reconstructed  by  the  numerical  algorithm  we  have  designed. 

I  gave  a  one  hour  talk  about  the  work  on  crack  detection  at  the  CBMS-NSF  conference 
on  Nondestructive  Evaluation  and  Inverse  Problems,  in  Lexington  in  June,  1995,  as  well 
as  a  plenary  talk  at  the  IFIP  conference  on  Modelling  and  Optimization  of  Distributed 
Parameter  Systems  with  Applications  to  Engineering,  in  Warsaw,  Poland,  in  July  1995.  In 
connection  with  my  talk  at  this  latter  conference  I  have  almost  completed  a  survey  article 
of  the  work  on  crack  detection  [2]. 

Jointly  with  a  graduate  student  (Donna  Fengya)  I  have  studied  the  effect  of  small 
inclusions  of  different  conductivity  inside  a  homogenous  conductor  [3].  We  have  derived 
an  asymptotic  representation  formula  for  the  voltage  potential  (including  a  convergence 
proof)  and  we  have  used  this  representation  formula  to  establish  Lipschitz  estimates  for 
the  locations  and  the  relative  sizes  of  the  inclusions.  On  one  hand  this  work  significantly 
extends  work  which  Avner  Friedman  and  I  carried  out  earlier  concerning  inclusions  with 
zero  or  infinite  condutivity,  on  the  other  hand  this  work  has  a  significant  computational 
component,  which  demonstrates  the  practical  viability  of  using  electrical  measurements 
to  find  such  small  inhomogenieties.  This  study  is  also  very  much  related  to  our  work 
to  demonstrate  the  viability  of  impedance  imaging  for  determining  corrosion  effects.  In 
this  context  we  have  made  good  progress  on  a  study  of  the  effective  current  boundary 
condition  for  a  metal  component  with  surface  corrosion  pits,  covered  by  a  dilutte  solution; 
this  study  has  involved  a  considerable  search  of  the  corrosion  engineering  literature.  A 
first  part  of  this  study  is  already  reported  in  a  paper  with  J.-M.  Xu  [4].  In  this  paper  we 
examine  the  questions  of  existence  and  uniqueness  of  solutions  for  a  model  (of  the  dilute 
solution)  based  on  the  socalled  Butler- Volmer  boundary  condition.  We  put  particular 
emphasis  on  the  regime  where  the  boundary  is  in  a  transition  to  passivity  (due  to  the  fact 
that  it  is  largely  covered  by^an  oxide  layer)  —  this  is  practically  and  mathematically  a  very 
interesting  regime.  It  is  worth  mentioning  that  we  have  developed  contacts  with  a  material 
science  group  at  the  EPFL  (headed  by  D.  Landolt).  We  hope  these  contacts  will  lead  to 
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experimental  work  concerning  the  transition  to  passivity.  In  a  second,  forthcoming  paper 
[5]  we  have  by  means  of  a  multiple  scale  analysis  derived  effective  boundary  conditions 
for  the  metal  component,  based  on  the  asumption  that  the  dilute  solution  is  very  shallow. 
Using  this  effective  boundary  condition  we  have  also  carefully  examined  the  question  of 
identifiabily  of  the  corrosion  damage.  Part  of  the  study  of  corrosion  imaging  has  been 
related  to  plates  -  we  have  for  instance  developed  (and  studied  the  performance)  of  a 
numerical  algorithm  to  image  corrosion  damage  on  one  side  of  a  plate  (the  interior)  from 
measurements  of  voltage  and  current  data  on  the  opposite  side  [6].  For  a  mathematically 
vaguely  related  problem  we  have  developed  a  method  to  reconstruct  the  interfase  between 
the  electrolytic  bath  and  the  aluminum  in  a  socalled  Hall-Heroult  cell  (a  standard  tool  to 
produce  aluminum).  This  cell  is  also  uthin  in  one  direction”  (it  is  very  shallow);  for  this 
problem  we  establish  an  estimate  of  the  error  between  the  interfase  reconstructed  from 
the  “shallow”  model  and  the  “true”  interface  (corresponding  to  the  higher  dimensional 
problem)  and  we  also  study  higher  order  approximations  [7].  The  electrostatic  data  is 
known  on  both  top  and  bottom  of  the  cell.  My  collaborator  on  this  project  (Michel 
Romerio  of  the  EPFL)  has  strong  ties  (as  a  consultant)  to  Alusuisse. 

The  book  project,  [8],  with  John  Sylvester  and  Gunther  Uhlmann  (entitled:  Electrical 
Impedance  Tomography)  has  progressed  steadily,  although  somewhat  slower  than  we  had 
initially  hoped.  One  of  the  reasons  for  the  delay  is  our  decision  to  include  a  section  with 
the  very  recent  results  about  two  dimensional  uniqueness  (due  to  A.  Nachman)  as  well  as 
a  section  (somewhat  longer  than  originally  envisaged)  dealing  with  the  most  recent  results 
about  crack  detection.  We  expect  to  have  a  complete  manuscript  finished  in  the  spring 
of  1996.  With  John  and  Gunther  I  have  also  organized  an  AMS  summer  conference  on 
Electrical  Impedance  Imaging,  which  took  place  in  Seattle  in  July  1995. 

The  work  on  first  order  corrections  to  homogenized  eigenvalues  has  continued  succes- 
fully.  We  study  the  possible  limits  of  (X€  —  X)/e  as  e  — >  0.  Here  A6  is  an  eigenvalue 
of  a  periodic,  composite  medium  with  material  parameters  a(x/e)  (e  is  the  microscopic 
length  scale)  and  A  is  the  corresponding  eigenvalue  of  the  effective  medium  with  constant 
material  parameters,  A.  We  call  these  limits  (A*)  first  order  corrections,  since  they  give 
rise  to  asymptotic  formulae  X€  —  A  +  eA*  +  o(e).  Jointly  with  a  graduate  student  (Shari 
Moskow)  I  have  succeded  in  proving  that  for  a  polygonal  domain  whose  sides  have  rational 
slopes  (relative  to  the  periodic  microstructure)  there  is  generally  not  just  one  correction 
but  frequently  a  continuum  of  possible  corrections  [9],  [10].  For  the  case  of  Dirichlet  bound¬ 
ary  conditions  this  result  had  already  been  established  during  the  last  reporting  period  - 
much  to  our  surprise  we  found  (during  this  reporting  period)  that  the  same  result  holds 
for  Neumann  boundary  conditions. 

We  provide  a  very  explicit  representation  formulas  for  the  set  of  possible  first  order 
corrections.  Whether  a  limijt  (of  ( X€n  —  X)/en)  is  attained  for  a  particular  sequence  en  — >  0, 
and  exactly  what  is  the  value  of  this  limit,  depends  on  the  interaction  between  the  pe¬ 
riodic  microstructure  and  the  macroscopic  boundary,  dCl.  Our  analysis,  which  is  quite 
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intricate,  involves  a  combination  of  error-estimates  for  approximation  of  eigenvalues,  a 
careful  study  of  the  minimal  smoothness  properties  required  to  derive  higher  order  esti¬ 
mates  of  homogenization  convergence  as  well  as  a  very  detailed  analysis  of  the  associated 
boundary  layers.  In  the  case  of  Neumann  boundary  conditions  we  also  make  extensive  use 
of  duality.  The  techniques  developed  have  immediate  implications  beyond  the  eigenvalue 
problem,  for  instance  for  the  boundary  layers  that  inevitably  appear  in  connection  with 
many  homogenization  problems.  It  is  exactly  because  of  this  close  link  to  boundary  layers 
that  we  were  initially  inclined  to  think  that  Neumann  boundary  conditions  would  lead 
to  behavior  significantly  different  from  that  exhibited  with  Dirichlet  boundary  conditions 
(they  do  in  one  dimension  -  but  not  in  two  dimensions!)  Some  of  results  we  have  proven 
were  conjectured  in  earlier  work  with  Fadil  Santosa  .  We  are  currently  in  the  process  of 
carrying  out  a  number  of  two  dimensional  calculations  of  the  first  order  eigenvalue  correc¬ 
tions  and  the  true  eigenvalues  (with  oscillatory  material  parameters).  This  work  (which 
will  be  contained  in  [11])  represents  a  significant  extension  of  our  earlier  calculations. 

My  investigation  concerning  certain  elliptic  inverse  problems,  originating  from  Toka- 
mak  fusion,  has  continued.  In  its  simplest  form  this  investigation  concerns  the  two  di¬ 
mensional  Laplace  operator.  The  problem  is  then  to  identify  /  >  0  of  the  equation 
A u  +  f(u )  =  0  in  Cl,  from  knowledge  of  the  boundary  flux  ^ u ,  corresponding  to  homo¬ 
geneous  Dirichlet  data.  Extending  our  previous  work,  Elena  Beretta  and  I  have  succeeded 
in  proving  that  if  the  domain  has  a  proper  corner  (with  angle  different  from  7r)  and  if  it 
is  apriori  known  that  /( 0)  =  0,  then  it  is  possible  to  recover  all  the  derivatives  of  /  at 
0  from  knowledge  of  -§^u\au,  [12].  If  it  is  not  known  that  /( 0)  =  0  then  the  same  result 
applies,  except  we  have  to  exclude  a  possibly  countable  number  of  angles  6  <  ir/2.  We 
believe  this  latter  exclusion  is  for  technical  reasons  alone,  and  that  the  recovery  should  still 
be  possible  in  the  case  of  these  angles.  It  is  not  that  difficult  to  see  that  recovery  of  /( 0) 
and  /'( 0)  is  always  possible  if  there  is  a  proper  corner.  Indeed  a  relatively  straightforward 
application  of  the  hodograph  transformation  yields  that  if  recovery  of  /( 0)  is  impossible 
(and  /  is  C°°)  then  the  domain  must  have  a  G°°  boundary.  We  suspect  that  if  recovery 
of  any  of  the  derivaties  of  /  at  0  is  impossible,  then  the  domain  must  indeed  be  a  disk, 
but  we  have  so  far  been  unable  to  prove  this  in  any  generality.  I  have  recently  been  able 
to  establish  some  results  to  this  effect  for  affine  /.  It  is  well  known  that  corners  (of  angle 
7r/2)  occur  in  plasmas  of  Tokamaks  based  on  the  socalled  divertor  technology.  We  are 
currently,  using  methods  of  hodograph  transformation,  in  the  process  of  examining  the 
reconstruction  of  /  for  domains  that  are  perturbations  of  disks  and  we  have  obtained  some 
encouraging  partial  results.  I  gave  a  plenary  talk  about  this  work  at  the  annual  meeting 
of  the  SIAM  Activities  Group  for  Control  and  Systems  Theory,  in  St  Louis  in  April  1995. 
A  paper,  surveying  this  work,  is  also  in  preparation  [13]. 

During  the  reporting  period  I  have  given  invited  talks  at  several  universities,  for 
instance  at  Texas  A&  M,  University  of  Texas  Austin,  Cornell  University,  University  of 
Maryland  Baltimore  County,  University  of  Maryland  College  Park,  Ecole  Polytechnique 


3 


Federale  de  Lausanne,  Aalborg  University,  Aarhus  University,  Universitat  Gottingen  and 
Eidgenossische  Technische  Hochschule  Zurich. 
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